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ABSTRACT 



^ Aims. We test our capability of deriving stellar physical parameters of giant stars by analysing 
a sample of field stars and the well studied open cluster IC 4651 with different spectroscopic 
methods. 

Methods. The use of a technique based on line-depth ratios (LDRs) allows us to determine 
with high precision the effective temperature of the stars and to compare the results with those 
obtained with a classical LTE abundance analysis. 

Results, (i) For the field stars we find that the temperatures derived by means of the LDR 
method are in excellent agreement with those found by the spectral synthesis. This result is 
extremely encouraging because it shows that spectra can be used to firmly derive population 
characteristics (e.g., mass and age) of the observed stars, (ii) For the IC 4651 stars we use the 
determined effective temperature to derive the following results, a) The reddening E{B — V) of 
the cluster is 0.12 ±0.02, largely independent of the color-temperature calibration used, b) The 
age of the cluster is 1.2±0.2 Gyr. c) The typical mass of the analysed giant stars is 2.0±0.2Mq . 
Moreover, we find a systematic difference of about 0.2 dex in logy between spectroscopic and 
evolutionary values. 

Conclusions. We conclude that, in spite of known limitations, a classical spectroscopic analysis of 
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giant stars may indeed result in very reliable stellar parameters. We caution that the quality of 
the agreement, on the other hand, depends on the details of the adopted spectroscopic analysis. 

Key words. Stars: late-type - Galaxy: open cluster and associations: individual: IC 4651 - 
Techniques: spectroscopic 



1. Introduction 

The detailed study of stellar populations in our own Galaxy and in its neighborhoods has 
received a major impulse in the last years, thanks to the use of large telescopes coupled 
to multi-object high-resolution spectrographs. 

The spectroscopic determination of chemical abundances in the atmospheres of stars 
may greatly contribute to our knowledge of galaxies. In fact, once a set of chemical 
abundances is complemented with a stellar age, it is possible to assess the age-metallicity 
relation, and to invert the observed color-magnitude diagr am (CMP) th us obtaining the 



Star Formation History (SFH) of t he Galaxy (see, e.g., iTolstovl l2005l and the Large 



Magellanic Cloud case illustrated by 



Cole et al 



20051 ). When individual stellar ages are 



unknown, this analysis becomes more uncertain because of the so-called age-metallicity 
degeneracy, i.e. the fact that old metal-poor stars can occupy the same region of the 
CMD as young metal-rich objects. 

In this context, open clusters provide fundamental tools because each cluster repre- 
sents an homogeneous sample of stars having the same age and chemical composition. 
Moreover, they are very suitable for the investigation of several issues related to stellar 
and Galactic formation and evolution. In particular, young open clusters provide infor- 
mation about present-day star formation processes and are key objects for clarifying 
questions on galactic structure, while old and intermediate-age open clusters play an 
important role in linking the theories of stellar and galactic evolution. 

The main classical tool to study cluster properties is the color-magnitude diagram, 
which suffers from several uncertainties and intrinsic biases, such as the limited knowledge 
of the chemical composition of the stars and the degeneracy in deriving the distance 
(and therefore age) and reddening. As a consequence, the photometric analysis of open 
clusters alone might not be so conclusive for an accurate determina tion of ages, distances . 



metallicities, masses, color excesses, and temperatures, as shown bv lRandich et al.l (j2005l ). 
Thus, spectroscopic methods to determine the effective temperature of cluster members 
are efficient techniques that are independent of the cluster reddening. This reddening can 
then be obtained by comparing the spectroscopic results to the photometric ones. 



Send offprint requests to: K. Biazzo, e-mail: katia.biazzo@oact.inaf.it 
* Based on observations collected at the ESO telescopes at the Parana! and La Silla 

Observatories, Chile. 
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Effective temperatures can be determined by imposing the condition that the derived 
abundance for one chemical element with many lines in the spectrum (typically Fe i) 
does not depend on the excitation potentials of the lines (hereafter this technique will 
be called as "spectral synthesis" and the "spectroscopic temperature" derived in this 
way will be indicated TJ^^*-"). Another spectroscopic method is based on the ratio of the 
depths of two lines having different sensitivity to effective temperature. This line-depth 
ratios (LDRs) technique provides an excellent measure of stellar temperature (hereafter 



( Gray fc Johanson 



1991 



Strassmeier & Schordan 


2000 


Grav & Brown 


2001) 



In this paper we apply the LDR method to derive effective temperature of 
nearby evolved field stars with good Hipparcos parallaxes (accuracy better than 10%, 



da Silva et al 



20061) and of giant stars of the intermediate- age open cluster IC 4651. For 



both groups, the temperature was previously derive d photometrically by means of color 



2004; 



da Silva et al 



indic es and spectroscopically by abundance studies (|Pasauini et al 
2006). In this work, we compare T^g'^ with the previous photometric and spectroscopic 



determinations in order to check temperatures obtained with different methods, and to 
understand the temperature range in which each method can be used. This will verify 
how accurately physical parameters can be obtained by inverting the spectroscopic re- 
sults. Moreover, for the stars belonging to the open cluster IC 4651, we are able to derive 
in a "spectroscopic way" color excess, mass and age. This is an important test because, 
to our knowledge, this is the first time that these parameters are determined in a cluster 
by means of the line-depth ratio method. Moreover, the possibility of comparing our 
spectroscopic results on IC 4651 with those obtained by the classical fitting of the main 
sequence provides a unique opportunity to cross-check our inversion method. 



2. Sample selections and observations 



We have considered seventy-one evolved field stars and six giant stars belonging to the 
intermediate-age open cluster IC 4651. 

The sample of field stars has been selected and analysed by|l 



Setiawan et al 



(120041) to 



derive the radial velocity (RV) variations along the Red Giant Branch and to investigate 
the nature of the radial v e locity long-term variations. Subsequently, the same sample was 



studied by 



da Silva et al 



(j2006l ) for the determination of radii, temperatures, masses, and 
chemical composition with the aim to understand how the RV variability is related to 
the stellar physical parameters. 



The stars in the cluster IC 4651 have been taken from the sample of 



Pasquini et al 



(|200J). The evolutionary status of the analysed stars is indicated by their position in the 
CMD shown in Fig. [1] in which our stellar sample is marked with arrows. 
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Fig . 1. CMP of IC 4651 . V and v — y are taken from 
and 



Meibom et al 



Pasquini et al 



Anthony- Twarog et al 



( 20021 ). respectively. The stars analysed by 



Meibom et al 



mm 



2002) and 



2004 ) are indicated by different symbols, while our ones are marked with 



arrows. 



The details of the obser vations, data reductio n, and instrumentation are described in 
( 2004 ) and iPasquini et al.l ([200J) . Their main characteristics are briefly 



Setiawan et al 



summarized here. T he spectra of the fie ld stars were acquired with the FEROS spectro- 
graph {R = 48 000. Ixaufer et al.lll999l ) at the ESO 1.5m-telescope in La Silla (Chile), 



while the IC 4651 spectra were acquired with the UVES spectrograph {R ~ 100 000, 



Dekker et al 



20001 ) at the ESO VLT Kueyen 8.2m-telescope in Cerro Paranal (Chile). 



The FEROS spectrograph has a wavelength coverage between 3700 and 9200 A, while 
the spectral region of the U VES observations cov ers the 5800-6800 A range (we have 



selected the data acquired bv lPasquini et al 



20041 with the CD#3 since they covered the 



spectral region of our interest). The signal-to- noise [S/N] ratio in all cases is higher than 
150/pixel and 160/pixel for the FEROS and UVES spectra, respectively. 



3. Data analysis 

3.1. Effective temperature 

The spectral regions covered by FEROS and UVES contain a series of weak metal lines 
which can be used for temperature determination with the LDR method. Lines from 
similar elements such as iron, vanadium, titanium, but with different excitation potentials 
(x) have different sensitivity to temperature. This is because the line strength is a function 
of temperature and, to a lesser extent, of the electron pressure (and consequently, the 
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surface gravity). This sensitivity arises from the exponential and power dependence on 
temperature in the excitation and ionization processes. Different lines are formed at 
different depths of the stellar atmosphere and therefore one would expect that lines from 
levels with high excitation potentials are formed deeper in the atmosphere and hence at 
higher temperatures. For this reason, it is better to choose pairs of lines with the largest 
X-difference. Between 6180 A and 6280 A there are several l ines of this type whose 



depth ratios have been exploited for temperature calibrations (jGrav fc Johanson 



Grav &: Brown 



2001 



Catalano et al 



2002 



Biazzo et al 



1991 



2007a!), to study t he rotational 



modulation of the average effective temperature of magnetically active stars (IFrasca et al 



2005 



Biazzo et al 



with phases 



2007b) or to investigate in Cepheid sta rs the pulsational variations 



Kovtvukh fc Gorlova 



weak line pairs for which 



200C; 



Biazzo et al 



Biazzo et al 



2004) . In particular, we chose 15 



(j2007al ) made suitable calibrations. These lines 



are of elements belonging to the iron group that have different temperature sensitivity. 
Furthermore, the S/N ratio around 150 for all the spectra makes them very suitable 



degree is expected 


Grav & Johanson 


1991; 


Grav Sz Brown 


2001; 


Kovtvukh et al. 


2006 


Biazzo et al. 


2007a 


. 



3.1.1. Sample of field stars 

All the field stars analysed in this paper have low rotational velocity (v sin i < 7 km 
s^^). The spectral resolution of FEROS is quite similar t o that of the spectro graph 
ELODIE {Observatoire de Haute- Provence, France), used by 



Biazzo 



Biazzo et al 



et al. 



( 2007a[ l. As a 



2007al ) at wsini < 5 



consequence, the LDR— Tcff calibrations developed by 
km have been applied to derive the effective temperature {T^g'^). Our star sample 
is comprised of sixty-seven giants and four stars with log q > 4.0 (n amely HD2151, 



da Silva et al 



HD16417, HD26923 and HD62644), as derived by 
the LDR— Tcff calibration for main sequence stars obtained by 
been used. 

We would like to note that 



( 2006 ). For the latter 



Biazzo et al 



(|2007al ) has 



Biazzo et al 



(2007a[) adopted the color-temperature rela- 



tionship obtained bv lGrayl (j2005l ). which does not include corrections due to metallicity. 



We shall come back to this point in the discussion of the Sections 14.1.11 and 14.2.11 
3.1.2. IC 4651 

For the giant stars of IC 4651, the UVES spectra have R = 100 000, significantly higher 
than the ELODIE {R — 42 000), so, in order to apply the calibrations, it was necessary 
to degrade the UVES spectra to the resolution of ELODIE. The analysed six giant stars 
have rotation al velocities v s' m i in t he range 0-15 km s~^, thus the LDR— Tcft calibrations 



developed by 



Biazzo et al 



(|2007a[ ) for the appropriate tisini have been applied to de- 
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rive the effective temperatures (T^^^). For the subgiant E95, an interpolation between 
the main sequence and the giant cahbrations was apphed. Consequently, the effective 
temperature uncertainty for this star should be higher than for the giants. 



3.2. Surface gravity 

We can compute the "photometric" surface gravity using the relation log 5 = log^Q + 
log(M/MQ) + 41og(roff/rofj 0) — log(L/LQ), where gQ is the solar surface gravity, while 
M, Tcff and L are the stellar mass, photometric effective temperature, and luminosity in 
the respective solar units. 



3.3. Reddening 



Seve ral determinations o f the photometric reddening of IC 4651 exist in literature (see. 



e-g-, 



Meibom et al 



2002I . and reference therein). Once the spectroscopic temperatures 



have been derived, we can estimate for each star the intrinsic color index {B — V)o by 
inverting the {B — V)o — Tcff relation and to compute the color excess E{B — V) from 
the observed and reddened {B — V) color. This can be obtained for each star separately, 
therefore the statistics on all stars can tell us not only the reddening of the cluster, but 
also provide a check of the goodness of the method applied. 



3.4. Mass and Age 



Mass and age are esti mated using a slightly modified ver sion of the Bayesian estimation 

( 20061 ) proposed 



da Silva et al 



method conceived bv I J0rgensen fc Lindegrenl (|2005r ). as 

(see references therein for more information). In short, given the stellar absolute mag- 
nitude, effective temperature, metallicity, and the associate errors, we can estimate the 
probability that such a star belongs to each small section of a theoretical stellar isochrone 
of a given age arid meta llicity. In particular, we have considered the isochrones developed 



by 



Girardi et al 



(j2000r ). Then, the probabilities are summed over the complete isochrone, 
and hence over all possible isochrones, by assuming a Gaussian probability of having the 
observed metallicity and its error, and a constant probability of having stars of all ages. 
The latter assumption is equivalent to assuming a constant star formation rate in the 
solar neighbourhood. In this way, at the end, we have the age probability distribution 
function (PDF) of each observed star. PDFs can also be obtained for any stellar property, 
such as initial mass, surface gravity, intrinsic colour, etcQ. 

Although a full discussion of the PDF method is beyond the scope of this paper, 
we note the following. The method, with just some small differences, has already been 



^ A Web version of this method 
http://stev.oapd.inaf.it/~lgirardi/cgi-bin/paraml 



available at the URL 
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tested on both main sequence stars (|Nordstrom et al 



2004) and on giants and subgiants 



(jda Silva et al.ll2006l ). Ages of dwarfs turn out to be largely undetermined by this method, 
due to their very slow evolution while on the main sequence. Ages of giants turn out to 
be well determined provided that the effective temperat ure and the parallax (absolute 
magnitude) are measured with enough accuracy. In fact. Ida Silva et al.l (|2006f ) find that 
stars with errors of 70 K in Tcff, and less than 10% errors in parallaxes, have ages 
determined with an accuracy of about 20%. These errors become larger for particular 
regions of the CMD, for instance on the red clump region, where stars of very different 
age and metallicity become tightly clumped together, and where in addition there is 
a superposition of red clump stars and first-ascent RGB ones. On the other hand, the 
best age determinations are expected on the subgiant branch region of the CMD, where 
evolutionary tracks of different masses separate very well from each other. To some extent, 
also the lower part of the red giant branch provides good age determinations, since in 
this CMD region there is no superposition with other evolutionary stages and stars above 
a given absolute magnitude can only be fitted by stars below a given mass (and above a 
given age). 



Of course, the above-mentioned errors do not include the systematic errors that may 
hide in the stellar evolutionary tracks and isochrones used. In fact, in any set of evolution- 
ary tracks the position of the giants depends on the choice of mixing length parameter 
(usually calibrated on the Sun), and to a much smaller extent on details of solving the 
atmosphere structure and on the interpolation of opacity tables. The only way to avoid 
this kind of uncertainty would be the building of tracks (and isochrones) in which the 
theory of energy transport and atmospheric structure are accurately calibrated on obser- 
vations of star clusters and binaries. This approach however is still too far to reach. Using 
the PDF method with different sources of evolutionary tracks could provide hints on the 
magnitude of systematic errors, but would not solve the situation because different sets 
of tracks in the literature do share the same assumptions and input physics. 



Therefore, at present we cannot reliably evaluate the systematic errors present in 
the PDF method. A simple check with two Hyadcs giants, where the PDF ages can be 
co mpared with the turn -off ones, hints to an accuracy of about 20% for the objects studied 



by 



da Silva et al 



( 20061 ). It would be extremely interesting to perform similar tests using 
giants in clusters of a wide range of ages, although we know that their distances (and 
turn-off ages) would be significantly more uncertain than the Hyades ones. 
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4. Results 



4.1. Sample of field stars 
4.1.1. Effective temperature 

In Table[T]we list names, Tjfp^ together with the photometric (Tj,^^*-*"^) an d spectroscopic 
(rpd'^ ^) temperatures, and as well as the stellar abundances retrieved by 



da Silva et al 



(|2006 ). The photometric temperature s were determined using the (B — V) — T^s rela- 



tionships of 



Alonso et al 



(1996 



1999f l for dwarf and giant stars. 



The comparison between these three temperatures is emphasized in Fig. [2] and dis- 
cussed here in the following. 

i. The temperatures derived by spectral synthesis (TJ^^^*-") and by LDR method (T^^^) 
are higher than the photometric ones by about 60-70 K, on average. This indicates 
that the photometric techniques tend to underestimate the temperature compared to 
the spectroscopic ones. This result seems to be more evident for the hottest stars in 
our sample. 

ii. The agreement between T^g^'^'^ an d the Tpd'^^ is bette r in the temperature range 



4000-5300 K, as aheady found by 



da Silva et al. 



(|2006f ) and displayed in the top- 



right panel of Fig. [2l The good agreement between T'™°'^ and TJj"^^ or T^fP^ in 
the 4000-5300 K temperature range is likely due to the fact that all these field stars 
are nearby and, consequently, the reddening, which is the parameter that mainly 
influences the photometric effective temperature determination, is negligible. The 
temperatures derived by spectroscopic methods are instead reddening- free. At higher 
temperatures, the behaviour of T^^^^'-'^ versus T^fjj'^^ or T^g*^ reflects the shape of 
the photometric calibration used for the temperature determination, 
iii. The temperatures obtained by means of the LDR method are in very good agree ment 



da Silva et al 



(20061), the 



in all the range with the spectroscopic values computed by 
average difference < T^p^ - TSPEC > ^^^^^^ ^bout 15 K with rms = 78 K. This rms 



can be considered as an estimate of the absolute error on the star temperature. 
IV. The TPHot 

versus T^fP^ plot shows larger scatter compared to T™'~'^ versus 



SPEC 



If we consider the whole temperature range, we obtain < T^p^ — T^fp'~''^ >= 73 K 
with the rms of 92 K. The mean difference < T^i^'^ - T™'^^ > is 58 K and the rms 
is 72 K. The fact that T™°t ^ldr gj^^^g higher scatter compared to T^fP°^ 
vs. rjf^^ can be due to the slight dependence of on metallici ty, especially 



Biazzo et al 



for star s with [Fe/H] significantly different than zero, as pointed out by 
(|2007a[) . In fact, if we plot the difference T^g^-T^g^'^ as a function of the metallicity 
(Fig. [3]), a residual dependence of the effective temperature on the metallicity emerges, 
mainly due to the metal-poor stars. If we disregard the five metal-deficient stars with 
[Fe/II]< —0.4 (asterisks in Fig. [2]), the average difference < T^^^-T^^^'^ > becomes 
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very low (5 K) with a value of 68 K as rms. If we instead take the linear relationship 
(full line in Fig. [3]) and correct T^fp^ for the residual dependence on the metallicity, 
the final residuals between the two temperatures T^g'^ and T^fj^^*^ is only 48 K. The 
agreement is remarkably good considering that in this small scatter are included all the 
possible differences, such as errors in the equivalent width and in the line-depth ratio 
measurements, errors in the LDR-Tcff calibrations and errors in the Tcs scale of the 
adopted standard stars. We think that the residual dependence of < T^fP^ — T^^^'^ > 
on the metallicity is due to the effect of [Fe/H] on LDRs that is a "second order" 
one compared to temperature and gravity, but can be neglected only for stars with a 
near solar metallicity. In order to adequately take into account this effect, a proper 
LDR— Tcff calibration should be done for a large sample containing many stars with 
well determined metallicity, spanning a wide range. 

One could object that the fact the two spectroscopic temperatures {TpP^ and T'^^g'^^) 
are in very good agreement should be expected, since the methods are very similar. 
However, there are a number of important differences that one should not, a priori, 
expect such a good agreement. 

— The two techniques are independent. One method simply measures several LDRs and 
takes advantage of the LDR~temperature calibrations of standard stars. These cali- 
brations use a {B — V) — Ten polynomial relationship to determine the temperature 
of the standard stars. The other method computes the temperature from the equiva- 
lent widths eliminating any dependence of the Fc i abundance on the line excitation 
potential and using appropriate atmospheric models. In order to compute the Fe i 
abundance, the process involves the use of stellar atmospheric models, the use of line 
strength [log gf), the determination of microturbulence (f) and gravity, a process 
which is not present in the LDR method. Only when all these steps are properly 
made, the results are similar. 

— We know that the lines are formed at different stellar atmospheric levels, but if one 
chooses pairs of weak lines of similar elements (such as elements of the iron group) , 
the abundance dependence is practically eliminated when the ratio between the two 
lines is considered. In the measure of the equivalent width, the abundance dependence 
is instead always present. 

— The effects of macroturbulence cancel out, to the first order, when the LDR is com- 
puted because they affect all lines in the same way. These effects are not present in 
the equivalent width measurements. 

— The continuum choice influences in a strong way the equivalent width determination. 
As a first approximation, this does not happen for the LDRs, because the continuum 
is practically the same for each line pair which are very close in wavelength. 
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— The rotational velocity shows its effect in the line-depth ratio computation, but can 
be taken into accou nt using appropriate calibrations at proper rotational velocity 



( Biazzo et al 



2007al) . The equival ent width measurement is not instead affected by 



rotational broadening (lGravll2005l ). 

Microturbulence is negligible when the depth ratio between a line pair is considered. 
The effective temperatures derived by spectral synthesis are instead obtained comput- 
ing the contribution of microturbulence. As a consequence, wrong microturbulence 
values (also of 0.1-0.2 km s^^) leads to wrong temperature determinations. 



Given all these differences, there are no a priori reasons why the effective temperatures 
resulting from the spectral synthesis and LDR methods had to be similar. 



4.2. IC 4651 



4.2.1. Effective temperature 

The star name and effective temperature as derive d with the LDR (Tj^^) and spectral 
synthesis (T^^f^^'-') methods (|Pasquini et al.l 12004) are listed in the first three columns 



of Table [2] Th e stars with prefi x "E" are from 



"MEl" is from 



Meibom 



Eggen 



( 19711 ). while the star with prefix 



(120001) . 

is always lower than T^^ff^*-^, by an amount between 70 and 90 K for most 
of the stars, a difference which is slightly larger than for t he field stars. One ca use is 



in the fact that IC 4651 is slightly metal rich ([Fe/H]~0.1, |Pasguini_et_al 



20041 ). This 



implies that we should apply a metallicity correction of 19 K to T^^^ to all the stars 
(cfr. previous Section). With this correction, the offset between the two temperatures 
becomes of 50-70 K, with the exception of E95. In the fourth column of Table [2] the 
values of the effective temperature {T^fP^"^) derived by the LDR method and corrected 
for metallicity are given. 

For the star E95 the difference is much larger, or 320 K. E95 is a subgiant star, as it 
appears from the CMD in Fi g.fl] This has two maj or effects in the present analysis. This 



is the evolved star for which 



Pasquini et al 



(200J) found the largest difference between 



photometric and spectroscopic temperature, and a full convergence of the solution could 
not be obtained. In addition, the LDR value quoted is a weighted average of the values we 
would obtain from the dwarfs calibration (^^5700 K) and the giants calibration ('^5250 
K). As previously mentioned, we do expect for this star a lower accuracy in its LDR 
temperature determination. 

The presence of only one hot subgiant in the sample does not allow us to derive 
firm conclusions, but it may indicate that before extending the method to stars in this 
portion of the HR diagram, further work is required. We note, however, that the change 
in temperature of 300 K for this star does not substantially change its derived age and 
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Table 1. Parameters of the field stars. 



Name 


(K) 


J.PHOT 

(K) 


'T^SPEC 
-'eff 

(K) 


[Fe/H] 


HD2114 


5259 


5180 


5288 


-0.03 


HD2151 


5973 


5785 


5964 


-0.03 


HD7672 


5121 


4957 


5096 


-0.33 


HD11977 


5024 


4901 


4975 


-0.21 


HD12438 


5236 


4888 


4975 


-0.61 


HD16417 


5872 


5729 


5936 


0.19 


HD18322 


4639 


4638 


4637 


-0.07 


HD18885 


4641 


4673 


4737 


0.10 


HD18907 


5143 


5056 


5091 


-0.61 


HD21120 


5203 


5026 


5180 


-0.12 


HD22663 


4514 


4792 


4624 


0.11 


HD23319 


4434 


4526 


4522 


0.24 


HD23940 


4962 


4782 


4884 


-0.35 


HD26923 


5985 


6126 


6207 


-0.06 


HD27256 


5077 


5012 


5196 


0.07 


HD27371 


4884 


4914 


5030 


0.13 


HD27697 


4901 


4876 


4951 


0.06 


HD32887 


4141 


4046 


4131 


-0.09 


HD34642 


4775 


4838 


4870 


-0.04 


HD36189 


5066 


4875 


5081 


-0.02 


HD36848 


4325 


4460 


4460 


0.21 


HD47205 


4604 


4777 


4744 


0.18 


HD47536 


4513 


4379 


4352 


-0.68 


HD50778 


4117 


4081 


4084 


-0.29 


HD61935 


4820 


4777 


4879 


-0.01 


HD62644 


5590 


5297 


5526 


0.12 


HD62902 


4168 


4230 


4311 


0.33 


HD63697 


4263 
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0.00 
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0.10 
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0.12 
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0.05 
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5025 
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0.07 
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0.13 
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5027 


4912 


4985 


-0.07 


HD99167 


4030 


3905 
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-0.36 


HD101321 
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4810 


4803 


-0.14 


HD107446 
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4145 


4148 


-0.10 


HDl 10014 


4414 


4429 


4445 


0.19 


HD111884 
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4271 


-0.06 
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5027 


4988 


5086 


0.09 


HDl 15478 


4252 


4293 


4250 


0.03 


HD122430 


4323 


4238 


4300 


-0.05 


HD124882 
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4293 


-0.24 


HD125560 


4395 


4443 


4472 


0.16 


HD131109 


4154 


4073 


4158 


-0.07 


HD136014 


4999 


4782 


4869 


-0.46 
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4694 


4654 


0.13 
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mass, since the evolutionary tracks in this portion of the HR diagram are substantially 
horizontal and the evolutionary phase is extremely fast. 

It is more difficult to prov ide a reliable estimat e for the systematic effects which 

(j2004l ). for instance, have in their sample 



could be present in our results. 



Carretta et al 



three stars (namely 27, 76, and 146 of their paper), with magnitudes and color indices 
very close to those of the stars E98 and E60 of our sample. In particular, their star 146 
corresponds to our E98. Yet, the spectroscopic temperatures t hey derive (4610 , 4620 and 
4720, respectively) are 180-290 K cooler than what found bv iPasauini et al.l (|2004l ) for 
E98 and E60, and 130-240 K cooler than derived by us with the LDR method. 

While a difference on a single star, of say 200 K, can be understood, a systematic 
effect of up to 200 K with respect to the LDR is more difficult to explain. T^g*^ for 
these stars have been obtained on high quality spectra and by using a large number 
of line-depth ratios (namely 15), which makes the measurements extremely robust. As 
far as the calibrations are concerned, our stars occupy the linear part of the LDR— Tcff 
relationships, far from saturation effects. We therefore consider that the T^fp^ measure- 
ment errors come mainly from the several LDR— T^^^^ relationships, which give different 
temperature values. Any systematic shift should therefore be dominated by the adop- 
tion of different temperature s cales. Since our calibrations are ultimately based on the 
{B — V")— temperature scale bv lGrayl ([2005|), we verified the systematic differences in the 
color range of our interest {B — V =0.80-1.13) for a number of calibrations present in 
literature. In Fig.[4l we can see that for giants of i? — ^ ~ 0.95 all the scales agree within 
100 K, and this agreement increases (~ 50 K) if we restrict ourselves to the most recent 
ones. We therefore tend to exclude s ystematic error s in Tp g larger than 100 K and argue 

(200J) as well as the reddening they 



that the temperature determined by 
have adopted are likely too low. 



Carretta et al 



4.2.2. Surface gravity 



To derive the photometric surface gravity of the giants studied in IC 4651 us i ng the 



Pasquini et al 



(200J) 



relation given in Section [3.21 the value of 1.8Mq suggested by 
as a lower limit of the mass of the turn-off stars was assumed, wh ile the photometric 
temperature T^^'~''^ was computed from the {h — y)~Tcs calibrations ( Alonso et al.lll999 ) 



Pasquini et al 



taking into account the color excess E{b — y) — 0.091 and [Fe/II]=0.10 
20041 ) . The luminosity L has been derived considerii ig the V abso lute magnitude of the 
stars, the bolometric correction [BC) tabulated by Flower ( 19961 ) as a function of Tcff 



and the solar bolometric magnitude M^^^y q = 4.75 (jCox 



2000f ). The absolute magnitude 



has been computed using the apparent visual distance m odulus V — My dete rmined by 

pOOOf ) with the 



the best fit of the is o chron es at Z=0.024 developed by 



Girardi et al 



CMD of 



Piatti et al. 



( 19981 ). For this isochrone fitting, the position of the red clump, 
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400 




Fig. 3. Metallicity dependence of T^^^ - T^i^'^. The five stars with [Fe/H]< -0.4 are 
showed with asterisks. Continuous line is the hnear fit to the points. 



as well as the main sequence and the turn-off, has been taken into account. This is an 
important point, because theoretical models predict that the absolute luminosity of the 
red clump stars depends fairly weakly on their chemical composition and age. Thus, by 
the fitting of the CMD we find V-My = 9.83 for an age of log(Age/yr) = 9.15±0.05 and 
E{B — V)=0.12. This corresponds to an absolute distance modulus of {V ~ Mv)o — 9.46. 
The appare nt distance m odul us we find just slightly excee ds the values of 9.7 and 9.8 



obtained by 



NissenI (|l988r i and lKjeldsen fc FrandsenI (jl99lh . respectively. 



In Table [2] the photometric gravit ies (log q^^^"*") are lis ted together with the spectro- 
scopic ones obtained (logg^^^^) bv iPasquini et al.l ( 20041 ) . whose gravities represent a 
lower limit because the Fei/Feii ionization equilibrium was not always reached. In Fig. [5] 
it is shown that the values of logg^^^*-" are on average « .2 dex higher than log g^^^"*" 



Carretta et al 



(200J) justify 



with differences also of up to 0.3 dex at lower temperatures, 
this effect, which they also found for their clump stars in IC 4651, as due to increased 
Fe II line blends a t lower temper a ture o r to variations of stellar atmospheric structure 



with temperature. 



da Silva et al 



(|2006l ) also found the same effect for their field stars, 
where the spectroscopic gravities are systematically overestimated by 0.2 dex in average, 
with a dependence on the effective temperature. 



4.2.3. Reddening 

T o deter i nine the redd ening of IC 4651, we have used the T^s — (B — V) calibrations 



of 



Cravl toom and of 



Alonso et 



([Fe/H]=0.10±0.03 for IC 4651, 



al 



1996 



Pasauini et al 



1999), wh ich includes the metallicity effect 



2004). Deriving the intrinsic color in- 
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Fig. 4. Examples of temperature calibrations published by several authors. 
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Fig. 5. Photometric gravities vs. spectroscopic gravities derived by 
for the giant stars studied in IC 4651. 



Pasquini et al 



( 20041 ) 



dex for all the stars from these relationships and comparing them with the observed 
colors, we obtain from the Gray's calibration E{B - V) ^ 0.111 ± 0.018, 0.164 ± 0.015 
and from the Alonso's calibration E{B - V) = 0.127 ± 0.016, 0.167 ± 0.017, in which 
the two values correspond to the results obtained with the two temperatures T^g'^ 
and TJj"^^, respectively. The color excesses obtained with the two calibrations are very 
similar and this improves if one considers the effective temperatures T^g'^'^ derived by 
the LDR method and corrected for metallicity effects. In this case the color excesses 
obtained with T^^P^'' is E{B - V) = 0.120 ± 0.016, which is near to the value ob- 
tained considering the Alonso's calibration. In addition, our color excess values are in 
good agreement with the resuhs of E{B - V) ^ 0.15, 0.10, 0.13, 0.10, 0.10 and 0.15 
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obtained by lEggenl (119711) 



(2002), 
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Nisse^(ll988D . iKieldsen fc Frandsenl (Il99lh . 



Anthony- Twarog fc Twaroe ( 20001) and 



Meibom et al 



2002), 



Pasquini et al 



Anthony-Twarog fc Twarog 



Meibom et al 



20041) ■ respectiyel y. Since 



2000D and 



Pasquini et al 



j 20041) 



compute E{b - y), the relation E{b - y) ^ 0.72E{B - V) (jCardeUi et all 119891 ) has 
been used. These results seem to confir i n our analysis of the Section 14.2.11 that the red- 

(200J) of E{B - y) = 0.083 ± 0.011 is likely too 



dening yalue deriyed by 
low. 



Garrett a et al 



4.2.4. Mass and age 



Since the PDF method allows us to deriye complete probability distribution functions 
separately for each analysed stellar parameter, we haye been able to deriye the mass of 
the six stars studied in IC 4651 and the age of the cluster. 

The results of the mass and age distributions based on the PDF method are plotted 
in Fig. [6] and Hsted in Tabled both for the temperatures T^P^ and T^i^'^- The mass 
and age distributions obtained for T^g'^'^ and T^^^^"^ are also shown for comparison. 
We find tha t the a yerage mass of the six giant stars studied in IC 4651 is 2.0±0.2 Mq. 



Pasquini et al 



(|2004l ) use stellar models of 1.8Mq to reproduce the lithium abundances 
as a function of temperature simultaneously for ten post-turnoff and giant stars (see their 
Fig. 5). If we instead consider only the six giant stars, the lithium abundance after the 
dredge-up is 0.2 ^ log N {Li) S 2.0. This lithium abundance can be explained by means 
of a 2.0Mq rotating model with initial rotational yclocity Vrot of 110 km s~^ and wit h 
a subsequent modest braking which leads to retain more lithium (jPalacios et al.l 12003( 1. 
This mass yalue is consistent with our determination. 

The cluster age found with the PDF method is 1.2±0.2 Gyr. This agrees with the 
yalue of log(Age/yr)=9.15±0.05 (1.40 ± 0.15 Gyr) that we determined yia isochrone 
fitting (Section 14.2. 2p , but is significantly lower tha n th e yalues of 2.4ifc0.3 Gyr and 
1.7±0.15 Gyr found by I Anthony-Twarog et all (|l988l ) and Meibom et all ([2002), respec- 
tiyely. Discussing the origin of these differences is beyond the scope of this paper. We just 
recall that the main differences between our and their isochrone fitting is that we consider 
also the position of the red clump in the GMD, additionally to the position of the main 
sequence and turn-off stars. The red clump proyides a strong constraint to the apparent 
distance modulus that can be used in the isochrone fitting, thus better constraining the 
cluster age via the turn-off. 

Clearly the errors associated to our ages and masses should be intended as internal 
errors. The fact of using different stars average out the errors on the single objects, 
however, all possible effects due, for instance, to a temperature scale shift and to the 
use of a particular set of evolutionary tracks are not considered. Systematic effects will 



change our results in a different way, depending on the position of the stars in the CMD. 
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Fig. 6. Mass and age distributions of the giant stars in IC 4651 for the different values 
of the effective temperature. 
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Table 2. Parameters of the giant stars analysed in IC 4651. 



Name 


rriLDR 
J off 


rpSPEC 
-l-cft 


rpLDRc 


My 


log/«°^ 


logs'™^ 


^LDR 


^SPEC 


Age^Di^ 


AgeSPE^c 




(K) 


(K) 


(K) 


(mag) 


(cm s^""") 


(cm s^""^) 


(Mq) 


(Mq) 


(Gyr) 


(Gyr) 


E12 


4926 


5000 


4945 


0.524 


2.50 


2.7 


2.464 


2.496 


0.698 


0.662 


E8 


4821 


4900 


4840 


0.867 


2.63 


2.7 


2.018 


2.156 


1.197 


0.999 


E60 


4833 


4900 


4852 


1.070 


2.69 


2.9 


2.013 


2.064 


1.175 


1.096 


E98 


4806 


4900 


4825 


1.080 


2.66 


3.0 


1.984 


2.061 


1.218 


1.098 


MEI 11218 


4906 


5000 


4925 


1.259 


2.75 


3.0 


2.027 


2.069 


1.115 


1.044 


E95 


5482 


5800 


5501 


2.153 


3.47 


3.5 


1.690 


1.582 


1.779 


2.166 



For instance, while the temperature is not of paramount importance in the Hertzsprung 
gap (star E95), it is the most critical parameter for Red Giant Branch stars. As far as 
the effective temperature scale is concerned, we tested the PDF method with the three 



stars 27, 76, and 146 studied by ICarretta et al.l |2004l ). with Toff = 4610,4620,4730 K, 



V = 10"^91, 10^^91, 10^94, respectively, and adopting their reddening {E{B-V) = 0.083) 
and their apparent distance modulus {V — My =10.15). The results are that for the 
stars 27 and 76 we obtain masses and ages well out of the accepted range (Age~ 3.4 
Gyr, M ~ 1.3M0), while for the star 146 the result (Age-- 1.7 Gyr, M ~ I.SA/q) is 
well compatible with our and within the range of the most accepted values available 
in literature. This shows that while a shift of ^^100 K doesn't change our analysis in a 
dramatic way, a difference of 200 K will produce unacceptable results. 
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5. Conclusions 

In this work we have derived accurate atmospheric parameters for field stars and for six 
giant stars in the open cluster IC 4651. 

For the former group of stars, with very low extinction and well determined physical 
pa rameters (Teg, l og q, £ ), we find a good agreement between the temperatures computed 



by 



da Silva et al 



( 20061 ) from spectral synthesis and those derived by us with the LDR 
technique. This gives strong support to the effectiveness of the LDR method in deriving 
effective temperatures. 

For the giant stars in the intermediate- age open cluster IC 4651, the LDR method 
allowed us to determine accurate effective temperature values, overcoming the problem 
of parameters degeneracy (in Tcff, logg, £,) encountered in the spectral synthesis method, 
and to derive the reddening of the cluster. We find tha t our value of E{B — V) is in 



agreement with previous works fe.g.. lPasauini et al 



derestimated by some other author (e.g. 
the stellar effective temperatures. 



2004 ), while it seems it has been un- 



Carretta et al .1120041) ■ who likely underestimated 



Thanks to the PDF method d eveloped bv lJorgensen fc Lindegren 
modified in 



da Silva et al 



20051 ) and slightly 



(|2006l ). we are able to find a mass of 2.0±0.2Mq for the giant 
stars and a cluster age of 1.2±0.2 Gyr, by adopting the effective temperature derived 
with the LDR technique. 

We conclude that our approach is well suitable to derive effective temperatures and 
reddening of evolved stars and clusters with a nearly solar-metallicity. The determination 
of very precise temperature is of great importance to derive stellar mass and age. For 
this reason, the LDR method could be inverted and used for stellar population studies 
in alternative and/or in addition to those based on photometric data. 
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